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建模策略 






























































从普通的股栗价格波动，到液体中悬浮微粒的位置所谓的布朗运动等分布领域广 
H. 简单的随机游动过程提供了一个很好的模型（至少在首次逼近的意义 上〉. 

































































































































余弦趋势 





























































































































































第 3 章 



最后，我们回到图表 1-1 中显示的洛杉矶年降兩量.在该序列中，我们没有发现序列依赖 
性的证据，但是现在寻找否定正态性的证据.图表 3-17 显示了该序列的正态分位数-分位数 
图.我们发现图形有明显的 弯曲. 一条穿过第一和第=•正态分位败的直线帮助指明图形对直线 
的偏离. 



本章的主要内容是描述、建模和估计时间序列的确定性趋势.最简单的确定性“趋势”是 
常数均值函数_估计常败均值的方法已经给出，但更重要的是研究了在不同条件下对估计量准 
确度的 评估. 然后，使用回归方法来估计时间上的线性或二次的 趋势. 接下来是关于周期性或 
季节性趋势建模的方法，还研究了所有这些回归方法的可靠性和有 效性. 最后一节开启了用于 










































































时间序列模型 




































平稳时间序列模型 

















图表 4-9 的散点图反映了极明 M 的负自相关. 






























































第 4 聿 


图表 4-19 显示了么 =1.5, 么= 一0.75下模拟的 AR (2) 序列的时间序 列图. 图表 4 -18中 
内 的周期性行为可以从具有相同周期（即360/30 = 12个时间 单位〉 序列的近似周期性变化清 
楚地反映出来.若0的单位是弧度，有时被称为 AR (2) 过程的准周期. 
_ ■ 表 4-19 », = 1.5, 6 = -0.75下 AB(2) 序荈的时 B 序 _ 



AR (2) 樓型的方差 

过程的方差 K 可以用 A 和 d 表示如下：对方程 （4.3.9) 两边求方差，得到 

r. = («+^)y.+2#,^y,+oJ (4.3.19) 

令灸=1,方程 （4.3.12) 给出了第二个关于 y 。 和 y, 的线性方程么 y ,， 结合方 
程（4.3.19〉得到 

y = (1-^)(1-~^4)-2^ = ( it ^) (1-^-^ (4 . 3 . 20 〉 


AR {2) 模型的少系数 

AR(2) 模型的一般线性过程表达式中的</>系数比 AR(1〉 的 复杂. 然而，我们可以用方 
程 （4.1.1) 写出 Y,, Y,-,, y, :的一般线性过程表达式，代人1=么1- 1 +先1-1+*,.如果 
我们令 A 的系数相同，得到递归公式： 

S 一^如=。 } (4. 3 . 21 ) 

— = 0 i ~ 2,3,—J 

通过递归公式，得到扣=1,咖=么， <h=^+*t， 以此类推.对于给定的 A 和通过 
这些关系式可以得到系数的数值解. 

还可以证明，在时，显式解为： 

. (4 - 3 ' 22) 
与前文一样， G 和(?：是 AR 特征方程的根的倒数.假如是复数根，方程 （4.3.22) 可以 
写作： 











































































































































E[log(y, 〉 ] 勿 log (/ <,), Var( log(Y,) 
对如下（泰勒） 展开式两边取期里和方差， 可得出上述 结果： 
log<Y,)»*log<^)+^^ 


柯.并且，如果该序列水平值大致呈指数变化， J 
接下来我们可能要进行一次差分，这将得出对取对数后的数据进行差分的 i 
百分比变动和对败 

，趋向于存在相对稳定的百分比变化. 


假设从一个时期到下一个时期， 

y , - (i+x,>y-, 

其中 100X, 是 y,-, 至 K 的百分比变化 （可 能为负数〉，则 

iog<y.) - iog(y^,) = log(^) = 

如果限制X,.例如， IXJC0.2 ( 


具体地说，假设 


,則可得 到-个 良好的通近， 


万千瓦小时为单位.较高的值与较低的值相比显示出更大的波动. 


n 表 5-8 典■月廑发电置 
































































































































































































































































































































































































































































































































































































































































































m(period.replace ) : qqline({ 


7.7 小结 

本章深人论述了 ARIMA 模抵的参数估计问埋，分别介绍了基于矩方法、各类最小二乘和极大 
化似然函数的估计准则，并且给出了各种估计貴的性质，还结合棋拟及真实时间序列数据加以例解. 
另外，还探讨并例解了应用自助法于 ARIMA 模型的问題. 

习题 

7.1 使用长度为100的序列，我们计算得到 n =0.8, r,=a5, r,=0.4, 7=2及样本方差 5. 假设 
—个带常数项的 AR(2) 棋型是适当的棋型，那么如何（简单）估计 A ，么， ft, of? 

7.2 假设下列败据得自某平稳过程，计算户，％,辟的矩估计： 6, 5, 4, 6, 4. 

7.3 如果 {Y,} 满足 AR(1) 模型，其中炎大约是 0.7, 那么要满足估计误差不超过 ±0.1, 在95% 
置信水平下需要多长的序列来估计 *=!Kl 

7.4 考虑己知均值是零的 ma < i 〉 过程.基于一个长度》=3的序列，观测到 y, -- o, y 2 = 

-1 ， y, = i/2. 

<a) 证明： （7 的条件最小二乘估计值是 1/2. 

(b) 求嗓声方差的估计.（提示，此简单情形无需使用迭代算法 .） 

7.5 对给定数据1=10, Y, = 9, y,=9.5, 拟合一个无常数項的 IMA(1, 1) 模型 • 

(a) 求 (? 的条件最小二乘估计值.（提 示： 先 做习® 7. 4.) 

(b) 估计 

















































































































型诊衡 








































































































































































































































































































































































第二项可写为： 


log/(y, 音 log2x -yijv ( --|-2 ^ 

其中； y(l|0>=0, v,-Var(y,). 


S 向量 z(0|0) 是；向镢，且 P<0|0> 是 z(o) — z(o|o)=z(o> 的 


CovC2 0 (i),i,0')]« i.i ■ 0，l“.，m — 1 









































































































































































i 时间序列回归模型 

p 列建模中纳人外部信息的想法.首先介绍模型包含着对时 
昏虑一类含有异常值影响作用的模型.相对于序列的正常行 

























































































































































































































































































































































































异方差时间序列模型 





















































































































































































































































































































































































































































































































































































































































































代败分 析 0 表明，周期图可视为一个“已平滑”的谱密度，所用谱窗为狄利克雷 （Dirichlet) 
核，见下式： 

注意，对所有形如/=_//» •的 傅里叶頻丰，都有 D(/) =0,因此该窗在这些频丰上不会有任何 
作用. 但是，图表 14-7 左边给出的 £>(/) 的图显示，主峰两侧有明显的“旁瓣”，这会导致非 
傅里叶頻串上的功率泄*进附近傅里叶频率上的功率中，如图表 14-6 所示. 

雉削是用于改善旁瓣问题的一种方法.锥削通过障低序列两端数据的 搞度， 使其数值慢慢 

趋于数据均 值零. 基本思路是，在计算有限长度序列的傅里叶变换时*降低边界效应.如果对 

序列进行锥削以后再来计算周期图，其效果就是对》»=100应用如图表 14-7 右边所示修正的狄 
利克雷核，如此旁瓣就被根本消除了. 


圈表 14-7 狄利克貢核与锥 M 后的狄《1克雷核 



頻率 _ mjL 


锥削最常用的形式基于余 弦钟. 用么代替原来的序列 y,， 

Y, = h,Y, (14. 5.3) 

其中的 A, 为余弦钟，例如由下式 给出： 

A,-|jl-cos[ 2at<t ； 0 ^ ) ]) <14.5.4) 

图表 14-8 左边给出了一个 n=100 的余弦钟的图 • 更常用的锥削由分裂余弦钟给出，即仅将余 
弦锥削应用于时间序列的极点而得到.分裂余弦钟锥削由下式 给出： 
fi-jl-cos^- 17 ^]), l<*<m 
A, = Jl , m + (14.5.^) 

n-m + l<t<« 

这称为一个 100p% 的余弦钟锥削，其中 p= 2 m/n. 图表 14-8 右边示出了一个10%分裂余弦钟 
锥削，其中 ”=100. 注意每端都有一个10%的锥削，总计有20%.实践中10%或20%的分裂 
余弦钟都很常用. 







































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































































附录 I 



















































































































































































































































































































































































